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ABSTRACT 

This  paper  considers  the  suitability  of  a  number  of  emerging  and  future  instruments  for  the  study  of  radio  recombination  lines  (RRLs) 
at  frequencies  below  200  MHz.  These  lines  are  of  interest  because  they  arise  only  in  low-density  regions  of  the  ionized  interstellar 
medium  and  because  they  may  represent  a  frequency-dependent  foreground  for  next-generation  experiments  trying  to  detect  Hi 
signals  from  the  Universe’s  Epoch  of  Reionization  and  Dark  Ages  (so-called  “21-cm  cosmology”  observations).  We  summarize 
existing  decametre-wavelength  observations  of  RRLs,  which  have  detected  only  carbon  RRLs.  We  then  show  that,  with  reference  to 
an  interferometric  array,  the  primary  instrumental  factor  limiting  detection  and  study  of  the  RRLs  is  the  areal  filling  factor  of  the 
array.  We  consider  the  first  station  of  the  Long  Wavelength  Array  (LWA-1),  the  LOw  Frequency  ARray  (LOFAR),  the  low-frequency 
component  of  the  Square  Kilometre  Array  (SKA-lo),  and  a  future  Lunar  Radio  Array  (LRA),  all  of  which  are  likely  to  operate  at 
decametre  wavelengths.  Key  advantages  that  many  of  these  arrays  offer  include  digital  signal  processing,  which  should  produce  more 
stable  and  better  defined  spectral  bandpasses;  larger  frequency  tuning  ranges;  and  better  angular  resolution  than  that  of  the  previous 
generation  of  instruments  that  have  been  used  in  the  past  for  RRL  observations.  Detecting  Galactic  carbon  RRLs,  with  optical  depths 
at  the  level  of  1(L3,  appears  feasible  for  all  of  these  arrays,  with  integration  times  ranging  from  a  few  hours  to  as  much  as  100  hr;  at 
optimal  frequencies  this  would  permit  a  Galactic  survey.  The  SKA-lo  and  LRA,  and  the  LWA-1  and  LOFAR  at  the  lowest  frequencies, 
should  have  a  high  enough  filling  factor  to  detect  lines  with  much  lower  optical  depths,  of  order  1CL4  in  a  few  hundred  hours.  The 
amount  of  RRL-hosting  gas  present  in  the  Galaxy  at  the  high  Galactic  latitudes  likely  to  be  targeted  in  Epoch  of  Reionization  and  Dark 
Ages  Hi  studies  is  currently  unknown.  If  present,  however,  the  spectral  fluctuations  from  RRLs  could  be  comparable  to  or  exceed  the 
anticipated  HI  signals. 
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1.  Introduction 

In  the  low-density  regions  of  the  Galactic  interstellar  medium 
(ISM),  free  electrons  can  be  captured  by  ions  at  very  high  quan¬ 
tum  numbers  (n  >  100).  As  the  atom  cascades  into  a  series  of 
successively  lower  ionization  states,  each  transition  produces  a 
radio  recombination  line  (RRL,  Gordon  &  Sorochenko,  2002). 
These  Rydberg  atoms  are  large,  often  having  macroscopic  di¬ 
mensions,  and  are  consequently  somewhat  fragile.  As  a  result, 
the  transitions  between  high  quantum  number  levels  are  a  sensi¬ 
tive  probe  of  the  environments  in  which  the  atoms  exist,  provid¬ 
ing  diagnostics  such  as  temperature,  density,  and  pressure,  and 
the  lines  can  be  used  to  constrain  the  size  of  the  regions  in  which 
the  atoms  occur.  RRLs  also  may  be  a  frequency-dependent  con¬ 
taminant  in  future  astrophysical  and  cosmological  observations 
using  the  highly  redshifted  Hi  21-cm  line  from  the  Epoch  of 
Reionization  (EoR)  or  earlier. 

Konovalenko  &  Sodin  (1980)  discovered  the  first  decametre 
RRL  at  a  frequency  of  26.13  MHz,  using  the  UTR-2  telescope 
to  observe  Cas  A.  Blake  et  al.  (1980)  subsequently  identified 
the  line  as  the  C63 1  a  line.  Since  then,  carbon  RRLs  have  been 
measured  towards  Cas  A  at  frequencies  from  14  to  1420  MHz 
(Payne  et  al.,  1994,  and  references  therein).  A  recent  study  by 


Stepkin  et  al.  (2007)  measured  a ,  (3,  y,  and  even  6  carbon-line 
transitions  at  frequencies  near  20  MHz.  The  lines  were  measured 
up  to  n  ~  1005,  which  is  the  highest  /; -bound  state  ever  detected. 

Observationally,  the  RRLs  resulting  from  high  Rydberg 
transitions  typically  occur  in  emission  at  frequencies 
above  200  MHz  (shortward  of  1.5  m)  and  have  line  widths 
dominated  by  Doppler  broadening.  Below  150  MHz,  however, 
the  excitation  temperature  of  the  atoms  approaches  the  typical 
gas  kinetic  temperature  and  the  lines  appear  in  absorption 
(Payne  et  al.,  1989).  In  between  is  a  poorly  studied  transition 
region,  in  which  the  lines  may  be  undetectable. 

There  are  a  number  of  new  telescopes,  either  under  construc¬ 
tion  or  in  the  design  phase  (Sect.  5),  that  will  be  operating  at  fre¬ 
quencies  below  about  200  MHz,  which  we  shall  henceforth  con¬ 
sider  “decametre  wavelengths.”  This  paper  considers  the  impact 
these  emerging  instruments  may  have  on  RRL  studies,  specif¬ 
ically  the  o'-transition  Carbon  RRLs.  We  review  aspects  of  de- 
cametric  RRLs  from  the  perspective  of  the  Galactic  interstellar 
medium  (Sect.  2)  and  cosmological  observations  (Sect.  3),  then 
in  Sect.  4  we  consider  the  detection  of  the  lines,  in  Sect.  5  we 
discuss  the  specifics  of  various  instruments  with  respect  to  RRL 
observations,  and,  in  Sect.  6  we  present  our  conclusions. 
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2.  RRLs  and  the  Galactic  ISM 

Pressure  and  radiation  broadening  as  well  as  Doppler  effects 
may  affect  the  RRL  absorption  line  profiles,  and  can  be  used  to 
derive  physical  properties  of  the  absorbing  cloud(s)  (Kantharia 
et  al.,  1998).  The  relative  strengths  of  the  a,/3,  and  y  lines,  which 
can  be  measured  simultaneously  at  different  quantum  numbers  n 
in  the  same  spectrum,  can  be  used  to  constrain  the  physical  prop¬ 
erties  of  the  absorbing  gas  further  (Erickson  et  al.,  1995).  While 
hydrogen,  helium,  and  carbon  lines  have  been  identified  in  emis¬ 
sion,  all  low-frequency  RRLs  known  to  date  are  high-/;  carbon 
transitions. 

Using  multi-frequency  model  fits  to  all  existing  data  on  the 
Cas  A  sightline  and  comparing  the  spatial  distribution  of  the 
higher  frequency  lines  to  that  of  Hi  and  molecular  clouds  on 
the  sightline,  Kantharia  et  al.  (1998)  concluded  that  these  RRLs 
most  likely  occur  in  diffuse  H  i  clouds  in  the  Perseus  arm,  which 
have  temperatures  of  order  75  K.  The  best  fit  physical  mod¬ 
els,  based  on  Cas  A  sightline  data  from  all  frequencies,  pre¬ 
dict  a  stronger  broadening  in  the  lowest  frequency  lines  than  has 
been  actually  observed.  According  to  Kantharia  et  al.  (1998), 
the  discrepancies  arise,  at  least  in  part,  from  difficulties  in  ac¬ 
curately  detecting  and  measuring  the  Lorentzian  wings  due  to 
the  weakness  of  the  lines  and  issues  with  accurate  spectral  base¬ 
line  removal.  Many  lines  reported  in  the  literature  have  been 
erroneously  fit  with  Gaussians  (e.g.,  Sorochenko  &  Smirnov, 
1990),  and  others  were  subjected  to  a  baseline  removal  pro¬ 
cess  which  significantly  diminished  the  fitted  line  strengths  and 
widths  (Payne  et  al.,  1994).  In  order  to  fully  realize  the  potential 
of  these  lines  in  measuring  physical  conditions  in  the  ISM,  sensi¬ 
tive  observations  with  robust  spectral  baselines  are  needed,  par¬ 
ticularly  at  the  lowest  frequencies.  With  their  very  high  signal- 
to-noise  observations,  Stepkin  et  al.  (2007)  argue  that  errors  in 
fitting  or  spectral  baseline  removal  alone  cannot  account  for 
the  line-width  discrepancy;  however  including  the  known  low- 
frequency  spectral  turnover  of  Cas  A  (e.g.,  Helmboldt  &  Kassim, 
2009)  lowers  the  amount  of  predicted  broadening  and  brings  the 
model  and  observations  into  agreement. 

Aside  from  the  Cas  A  sightline,  carbon  recombination  lines 
in  absorption  have  been  observed  on  a  variety  of  sightlines  in 
the  Galaxy  (Tables  1  and  2).  Most  observations  have  focussed 
on  sightlines  towards  known  bright  background  sources  or  pass¬ 
ing  through  gas-rich  regions  along  the  inner  Galactic  plane,  but 
the  detection  rates  have  been  low.  Only  two  broader  surveys 
for  the  lines  have  been  published.  Erickson  et  al.  (1995)  found 
carbon  RRL  absorption  complexes  on  roughly  30  sightlines  be¬ 
tween  Galactic  longitudes  of  340°  and  20°,  with  Galactic  lat¬ 
itudes  \b\  <  2°.  They  concluded  that  there  is  no  evidence  for 
pressure  broadening  in  their  data,  and  that  the  linewidths  re¬ 
flect  Doppler  broadening  and  the  presence  of  many  clouds  in 
their  4°  beam.  Kantharia  &  Anantharamaiah  (2001)  found  RRLs 
on  9  sightlines,  with  six  of  those  in  the  Galactic  longitude  range 
352°  <  £  <  17°.  Combining  their  data  with  that  of  Erickson  et 
al.  (1995)  and  information  from  observations  on  the  same  sight¬ 
lines  at  higher  frequencies,  they  conclude  that  the  clouds  likely 
have  sizes  of  2°— 4°  and  occur  in  photodissociation  regions. 

3.  RRLs  and  21-cm  Cosmology 

There  is  growing  interest  in  the  use  of  the  21-cm  hyperfine  tran¬ 
sition  of  H  i  as  a  cosmological  and  astrophysical  probe  of  the 
distant  Universe.  For  comprehensive  reviews,  see  Furlanetto  et 
al.  (2006)  and  Pritchard  &  Loeb  (2008).  For  a  portion  of  the 
Universe’s  history  (1100  <  z  <  7),  the  baryonic  content  of  the 


intergalactic  medium  (IGM)  is  dominated  by  Hi,  and  interac¬ 
tions  between  the  adiabatic  expansion  of  the  IGM  and  various 
heating  sources  (e.g.,  first  stars,  first  black  holes,  dark  matter 
decay)  can  produce  either  an  absorption  or  emission  Hi  signal 
relative  to  the  cosmic  microwave  background  (CMB). 

At  least  three  distinct  epochs  have  been  identified,  which 
have  corresponding  spectral  windows  for  the  redshifted  H  i  sig¬ 
nal. 

1.  The  Dark  Ages  -  100  <  z  <  30  (15  MHz  <  v  <  50  MHz). 
The  H  i  gas  is  expanding  and  cooling  at  a  faster  rate  than  the 
CMB,  and  collisions  within  the  gas  drive  the  kinetic  temper¬ 
ature  of  the  gas  below  the  CMB  temperature,  resulting  in  an 
absorption  feature.  This  epoch  is  thought  to  occur  before  the 
first  stars  form,  and,  absent  any  other  heating  sources  (e.g., 
energy  injection  by  decaying  dark  matter,  Furlanetto  et  al., 
2006),  the  evolution  of  this  signal  should  depend  only  on 
cosmological  parameters. 

2.  First  star  formation  -  30  <  z  <  15  (50  MHz  <  v  <  90  MHz). 
As  the  first  stars  form,  they  flood  the  Universe  with  Lyman-cr 
photons.  These  strongly  couple  the  kinetic  and  spin  temper¬ 
atures  of  Hi  and  produce  a  second,  deep  absorption  feature. 
Detection  of  this  signal  would  constrain  both  the  epoch  and 
luminosity  function  of  the  first  stars. 

3.  Epoch  of  Reionization  (EoR)  -  15  <  z  <  7  (90  MHz  <  v  < 
200  MHz).  As  stars  and  the  first  black  holes  heat  the  gas,  its 
temperature  rises  above  the  CMB  temperature.  The  absorp¬ 
tion  feature  turns  into  an  emission  feature,  which  persists  un¬ 
til  reionization  completes  and  the  signal  cuts  off.  Detection 
of  this  feature  would  constrain  the  duration  of  the  EoR  and 
the  luminosity  function(s)  of  ionizing  sources  present. 

Estimates  of  the  strength  of  these  H  i  signals  is  both  epoch  and 
model  dependent  but  are  in  the  range  10-100  mK. 

One  of  the  key  challenges  to  detecting  any  of  these  signals 
are  the  various  foregrounds.  These  foregrounds  include  radio 
frequency  interference  (RFI)  from  terrestrial  radio  transmitters, 
ionospheric  phase  corruptions,  the  Galactic  synchrotron  emis¬ 
sion,  and  emission  from  extragalactic  sources.  The  various  spec¬ 
tral  windows  for  these  H  i  signals  also  include  frequencies  where 
RRLs  have  been  observed  on  Galactic  plane  sightlines.  It  is  cur¬ 
rently  unknown  at  what  strength  RRLs  might  exist  on  sightlines 
at  the  high  Galactic  latitudes  where  most  of  the  cosmological 
studies  will  be  made. 

A  standard  assumption  for  the  removal  of  astrophysical  fore¬ 
grounds  (Galactic  synchrotron,  extragalactic  sources)  is  that 
their  spectra  are  smooth  (Harker  et  al.,  2009).  While  21-cm 
cosmological  observations  will  clearly  be  conducted  at  high 
Galactic  latitudes,  in  order  to  minimize  the  contribution  of  the 
Galactic  synchrotron  emission,  even  residual  effects  from  RRLs 
might  still  vitiate  the  measurements. 

As  an  illustration  of  the  potential  impact  of  RRLs  on  these 
experiments,  consider  a  set  of  RRLs  having  an  optical  depth 
r~5x  10~4,  characteristic  of  what  has  been  observed  to  date 
(Table  1).  We  assume  no  continuum  opacity,  as  is  reasonable 
for  observations  at  frequencies  v  >  10  MHz  at  high  Galactic 
latitudes.  The  continuum  temperature  Tc  will  be  the  combina¬ 
tion  of  the  Galactic  synchrotron  emission  and  the  extragalactic 
background  but  is  likely  to  be  dominated  by  the  former  (Bridle, 
1967).  The  temperature  in  the  line  is  then  (Tc  +  T)t,  for  a  tem¬ 
perature  T  of  the  gas  hosting  the  RRLs  and  employing  the  usual 
expansion  for  small  optical  depths.  We  now  consider  the  two 
cases,  Tc  »  T  and  Tc  ~  T. 

Using  150  MHz  as  a  characteristic  frequency  for  detection 
of  the  H  i  signal  from  the  EoR,  the  coldest  parts  of  the  sky  have 
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a  temperature  Tc  «  150  K.  In  the  Galactic  plane,  the  gas  hosting 
the  RRLs  has  T  ~  100  K  (Sect.  2),  and  if  there  is  such  RRL- 
hosting  gas  at  high  Galactic  latitudes,  then  clearly  T  ~  Tc.  In 
this  case,  spectral  fluctuations  between  RRLs  and  the  continuum 
of  order  (Tc  +  T)t/Tc  ~  2 r  should  be  expected.  Conversely,  if 
one  considers  75  MHz  as  a  characteristic  frequency  for  the  H  i 
signal  from  the  first  star  formation,  the  coldest  parts  of  the  sky 
have  a  temperature  Tc  «  1000.  Assuming,  as  before,  T  =  100  K, 
spectral  fluctuations  of  order  r  should  be  expected. 

The  magnitude  of  the  spectral  fluctuations  expected  from  the 
Hi  signal  varies  with  the  epoch  being  considered,  but  a  typical 
range  of  values  is  10  5— 10  4.  If  there  is  RRL-hosting  gas  at  high 
Galactic  latitudes,  the  spectral  fluctuations  arising  from  RRLs 
could  be  comparable  to  or  exceed  that  from  the  desired  H  i  sig¬ 
nals. 

4.  Detection  of  RRLs 

We  begin  by  reviewing  the  observational  requirements  for  de¬ 
tecting  RRLs.  From  the  radiometer  equation,  the  sensitivity  of 
an  experiment  designed  to  detect  RRLs  can  be  expressed  as 


and  therefore  different  angular  resolutions,  may  help  resolve  the 
extent  to  which  a  mismatch  has  been  an  issue  in  previous  obser¬ 
vations. 

The  system  temperature  at  low  radio  frequencies  is  typically 
dominated  by  the  synchrotron  emission  from  the  Galactic  plane 
itself  and  thus  varies  with  both  position  and  frequency.  As  can 
be  seen  in  Table  1,  most  of  the  Galactic  lines  are  very  weak, 
with  optical  depths  r  <  1 0  5 .  Even  a  5cr  detection  (i.e.,  m  =  5, 
equation  2)  therefore  implies  AT /Tsys  <  2  X  10  4.  Any  survey 
for  these  lines  should  be  targeted  to  reach  at  least  this  sensitivity 
level. 

The  lines  themselves  vary  in  velocity  width  from  10 
to  100  km  s^1,  suggesting  that  a  velocity  resolution  of  roughly 
2  km  s  1  (~  700  Hz  at  100  MHz)  is  appropriate  for  blind  detec¬ 
tion  of  the  narrowest  lines.  However,  even  observations  targeting 
the  broader  lines  will  benefit  from  high  spectral  resolution  to  aid 
in  the  identification  and  excision  of  RFI.  Because  the  lines  them¬ 
selves  are  quite  weak,  even  weak  RFI  must  be  carefully  removed 
in  these  observations. 
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5.  Prospects  for  Future  Observations 


where  AT  is  the  rms  noise,  Tsys  is  the  telescope  system  temper¬ 
ature,  /  is  the  filling  factor  of  the  telescope  or  array,  Armt  is  the 
integration  time,  Av  is  the  frequency  resolution,  /Vp()|  is  the  num¬ 
ber  of  polarizations,  and  Mines  is  the  number  of  simultaneous 
transitions  observed  that  can  be  folded  together. 

If  the  RRL  has  an  optical  depth  r,  then  an  “m- sigma”  detec¬ 
tion  implies  mAT  <  rTsys,  with  a  corresponding  integration  time 
of 

At-  =(  —  T"ys  f  1  (2) 
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We  have  introduced  the  beam  dilution  factor  (dilution  factor 
or  DF)  of  the  absorbing  signal,  which  is  relevant  for  the  case 
of  observations  of  lines  against  a  non-discrete  background  such 
as  the  Galactic  synchrotron  emission.  For  example,  estimates 
from  two  surveys  suggest  that  absorption  lines  seen  towards 
the  central  Galactic  plane  probably  arise  in  clouds  that  subtend 
between  2°  and  4°  on  the  sky,  but  are  unlikely  to  be  bigger 
(Kantharia  &  Anantharamaiah,  2001).  If  the  beam  of  the  antenna 
is  substantially  larger  than  4°,  it  is  likely  that  the  absorption  fea¬ 
ture  will  be  diluted  by  unabsorbed  background  radiation.  We  es¬ 
timate  DF  as  the  geometric  ratio  of  the  cloud  and  beam  sizes, 
and  DF  <  1.  The  dilution  factor  DF  enters  equation  (2)  as  a  mul¬ 
tiplier  of  AT,  so  A tint  K  DF-2.  This  factor  is  not  important  when 
the  background  source  is  discrete  and  thus  fully  covered  by  the 
absorbing  gas,  such  as  in  the  case  of  Cas  A  and  other  bright 
Galactic  or  extragalactic  sources. 

Most  previous  RRL  observations  have  been  conducted  at  rel¬ 
atively  low  angular  resolutions  (Table  2).  As  a  consequence,  a 
significant  uncertainty  has  been  the  extent  to  which  different 
source  regions  are  blended  within  the  telescope  beam.  One  of 
the  advantages  of  many  of  the  future  telescopes  (Sect.  5)  is  that 
they  will  have  better  angular  resolution  than  most  of  the  pre¬ 
vious  systems  used  for  RRL  observations. Thus,  we  include  the 
beam  dilution  factor  DF  as  a  generic  factor  to  account  for  a  po¬ 
tential  mismatch  between  the  angular  sizes  of  the  of  foreground 
absorber(s)  and  the  background  source(s)  in  the  beam.  Indeed, 
future  observations,  particularly  those  at  different  frequencies 


A  new  generation  of  low  frequency  instruments  have  either  the 
opportunity  to  contribute  to  Galactic  RRL  and  ISM  studies  or 
to  be  confused  by  the  lines  when  doing  other  studies  (Table  3). 
Improvements  in  both  sensitivity  and  resolution  will  help  further 
resolve  issues  of  physical  cloud  parameters  such  as  size,  pres¬ 
sure  and  temperature,  while  allowing  a  more  complete  census  of 
the  ionized  carbon  gas  in  the  Galaxy.  Moreover,  in  many  cases, 
one  of  the  prime  science  drivers  for  the  telescope  is  observations 
of  the  cosmological  H  i  signal,  so  that  careful  removal  of  fore¬ 
grounds,  including  RRLs,  is  paramount.  We  shall  consider  the 
performance  of  these  telescopes  in  turn. 

At  low  frequencies,  the  system  temperature  Tsys  should  be 
dominated  by  the  Galactic  synchrontron  emission,  and  all  of  the 
telescopes  will  have  the  capability  to  conduct  dual  polarization 
observations  (Mpoi  =  2).  Thus,  with  respect  to  the  performance 
of  these  future  telescopes,  most  notably  those  still  in  the  design 
phase,  equation  2  indicates  that  there  are  only  two  other  factors 
under  the  control  of  those  designing  the  telescope.  One  is  the 
filling  factor  /  of  the  telescope  and  the  other  is  the  number  of 
lines  Mimes,  which  is  coupled  to  the  total  bandwidth  that  the  sys¬ 
tem  can  handle. 

One  feature  common  to  all  of  these  future  systems  is  the  sig¬ 
nal  chain  will  be  largely  digital,  which  should  enable  significant 
improvements  in  the  stability  and  shape  of  the  spectral  band- 
passes.  Further,  all  of  the  telescopes  will  have  a  significantly 
enhanced  frequency  tuning  range  (with  typical  ranges  of  2:1 
or  better),  which  will  enable  many  more  lines  to  be  observed. 
Finally,  all  of  the  systems  that  we  will  discuss  will  be  dipole- 
based  phased  array  systems  in  which  the  dipoles  will  be  de¬ 
ployed  at  essentially  ground  level.  Such  systems  are  in  contrast 
to  single  dish  observations  in  which  the  feed  is  at  an  altitude  of 
tens  of  meters  or  more.  At  such  altitudes,  a  feed  is  susceptible 
to  many  low-level,  scattered  RFI  sources,  including  those  from 
over  the  horizon,  and  experience  has  shown  that  single  dish  RRL 
observations  are  rarely,  if  ever,  successful  during  the  day,  even 
on  days  when  the  general  RFI  level  should  be  quite  low  (e.g.,  the 
Easter  holiday).  Dipoles  deployed  at  ground  level  should  sup¬ 
press  RFI  sources  located  near  the  horizon  and  should  enable 
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more  sensitive  RRL  observations,1  a  hypothesis  generally  con¬ 
firmed  by  daytime  observations  with  the  dipole-based  UTR-2. 

5.1 .  Long  Wavelength  Array  (LWA) 

The  Long  Wavelength  Array  (LWA)  is  intended  to  be  a  next- 
generation,  high  angular  resolution  imaging  telescope  operating 
at  frequencies  between  20  and  80  MHz  (Ellingson  et  al.,  2009). 
Some  of  the  key  science  projects  are  the  high  redshift  Universe, 
including  searches  for  high-redshift  radio  galaxies  and  studies  of 
clusters  of  galaxies;  probing  particle  acceleration  in  supernova 
remnants,  radio  galaxies,  and  clusters  of  galaxies;  and  searches 
for  radio  transients.  The  telescope  would  be  composed  of  as 
many  as  50  dipole  phased  array  “stations,”  with  each  station  con¬ 
taining  256  dual  polarization  dipoles.  The  first  station  (LWA-1) 
is  currently  under  construction  and  is  planned  to  conduct  science 
observations  in  its  own  right. 

The  configuration  of  LWA-1  is  a  pseudo-random  distribution 
of  the  dipoles  (Kogan  &  Cohen,  2009),  aimed  at  minimizing  the 
sidelobes  from  the  phased-array  station.  Because  of  the  broad 
frequency  range,  the  filling  factor  of  the  station  is  frequency  de¬ 
pendent. 

Table  4  summarizes  our  estimates  for  the  required  integra¬ 
tion  times  to  detect  Galactic  plane  RRLs  with  LWA-1  over  its  op¬ 
erational  frequency  range.  We  assume  5<x  detections  of  roughly 
1 0  3  optical  depth  lines  at  roughly  2  km  s  1  velocity  resolution. 
We  have  also  assumed  an  absorber  size  of  4°,  which  is  compara¬ 
ble  to  the  station  resolution  at  the  higher  frequencies. 

This  first  station  of  the  planned  LWA  instrument  is  ideal  for 
RRL  detection  work  due  to  both  its  flexible  frequency  coverage 
and  its  high  filling  factor  at  lower  frequencies.  At  optimum  fre¬ 
quencies,  the  desired  sensitivity  may  be  reached  in  as  little  as 
6  hours. 

While  there  has  been  some  consideration  of  precursor  Dark 
Ages  observations  with  LWA-1,  the  prime  science  mission  for 
the  LWA  as  a  whole  does  not  include  cosmological  H  i  measure¬ 
ments.  Moreover,  it  is  not  clear  that  LWA-1  would  have  the  req¬ 
uisite  sensitivity  to  detect  the  Dark  Ages  signal.  Consequently, 
RRLs  are  not  likely  to  be  important  as  a  contaminant  in  any  other 
LWA-1  observations. 

5.2.  LOw  Frequency  ARray  (LOFAR) 

The  LOw  Frequency  ARray  (LOFAR)2  is  a  next-generation,  low 
radio  frequency  telescope  operating  at  frequencies  between  30 
and  240  MHz.  Its  frequency  coverage  is  split  between  a  high- 
band  antenna  (HBA,  110-240  MHz)  and  low-band  antenna 
(LBA,  [10]  30-90  MHz). 

One  of  its  key  science  projects  is  to  study  the  Epoch  of 
Reionization  using  the  high-band  antenna  at  redshifts  1 1.5  >  z> 
6.5  (115  MHz  <  v  <  190  MHz).  This  is  the  frequency  range  over 
which  the  RRLs  switch  from  absorption  to  emission;  at  some  fre¬ 
quencies  during  this  transition  they  may  be  undetected  (Payne  et 
al.,  1989). 

In  its  full  configuration,  LOFAR  has  a  very  low  areal  fill¬ 
ing  factor,  making  it  unsuitable  for  RRL  detection  and  study. 
However  using  the  “super-station”  core  of  6  stations  within  a 
300  m  diameter  area  or  using  a  single  station,  more  reasonable 
filling  factors  of  order  10%  can  be  achieved.  The  instrument 
has  flexible  frequency  coverage  and  it  is  possible  to  cover  large 

1  Presumably  the  suppression  of  RFI  should  be  less  important  for  the 
Lunar  Radio  Array  (Sect.  5.4). 

2  http://www.lofar.org 


bandwidths  (up  to  48  MHz)  simultaneously.  The  ability  to  fold 
together  many  lines  will  help  with  detection.  However,  the  fre¬ 
quency  setup  has  a  fixed  frequency  resolution,  corresponding  to 
about  5  km  s  1  velocity  resolution  at  the  lowest  frequencies;  thus 
the  narrower  RRLs  may  be  under-resolved. 

By  folding  together  hundreds  of  transitions  and  using  the 
super-station  core,  this  instrument  would  be  able  to  detect  10  3 
optical  depth  lines  in  a  few  hours  at  the  lowest  frequencies; 
longer  integrations  should  be  able  to  detect  these  and  fainter 
lines  at  higher  frequencies  or  with  less  folding  or  both.  Folding 
together  adjacent  RRLs  assumes  that  An  <sc  n,  and  thus  that 
the  transitions  are  essentially  identical.  This  assumption  is  un¬ 
likely  to  be  warranted  if  hundreds  of  lines  are  folded  together, 
and  the  resulting  line  profile  would  no  longer  be  useful  for  de¬ 
riving  physical  parameters  of  the  absorbing  gas.  The  planned 
EoR  experiment  will  integrate  deeply  enough  to  detect  RRLs,  if 
present,  at  the  higher  frequencies. 

5.3.  Square  Kilometre  Array  (SKA) 

The  Square  Kilometre  Array  (SKA)3  will  be  one  of  a  suite 
of  new,  large  astrophysics  facilities  for  the  21st  century  and 
will  probe  fundamental  physics,  the  origin  and  evolution  of  the 
Universe,  the  structure  of  the  Milky  Way  Galaxy,  and  the  forma¬ 
tion  and  distribution  of  planets.  Highlighting  the  SKA  Science 
Case  are  Key  Science  Projects  (KSPs,  Carilli  &  Rawlings,  2004; 
Gaensler,  2004),  one  of  which  is  “Probing  the  Dark  Ages  and 
the  Epoch  of  Reionization”  (Carilli  et  al.,  2004).  The  goal  is  to 
track  the  transition  of  the  originally  neutral  IGM  into  its  cur¬ 
rent  mostly  ionized  state  via  imaging  of  the  (highly-redshifted) 
Hi  line.  With  reference  to  the  three  21 -cm  cosmology  epochs 
(Sect.  3),  the  goal  for  the  SKA  is  to  cover  the  EoR,  and  as  far 
back  into  the  epoch  of  first  star  formation  as  possible  (20  <  z  < 
7,  70  MHz  <v<  200  MHz). 

The  configuration  of  the  low-frequency  component  of  the 
SKA,  the  so-called  SKA-lo,  remains  a  matter  of  active  study, 
but  the  current  specifications  forsee  approximately  25%  of  the 
total  collecting  area  within  a  diameter  of  1  km  and  approxi¬ 
mately  50%  of  the  total  collecting  area  within  a  diameter  of  5  km 
(Schilizzi  et  al.,  2007).  The  resulting  filling  factors  are  /  ~  0.5 
(/  ~  0.05)  within  the  central  1  km  (5  km). 

Scaling  from  the  results  in  Table  4,  the  central  portion  of  the 
SKA-lo  could  therefore  detect  RRLs  along  the  Galactic  plane 
in  integration  times  of  order  5  hr.  Given  the  primary  scientific 
driver  for  the  SKA-lo,  however,  it  is  unlikely  that  it  will  be 
conducting  deep  observations  along  the  Galactic  plane.  Rather, 
it  will  likely  conduct  its  deep  observations  (>  100  hr)  at  high 
Galactic  latitudes.  The  strength  and  distribution  of  RRLs  at  high 
Galactic  latitudes  is  currently  unstudied;  it  is  unknown  if  they 
will  be  present  at  all.  However,  even  if  the  RRLs  on  these  sight¬ 
lines  are  an  order  of  magnitude  weaker  than  currently  known 
lines,  with  r  ~  10~4,  equation  (2)  shows  that  they  would  still  be 
detected  in  integrations  of  order  100  hours. 

5.4.  Lunar  Radio  Array  (LRA) 

The  Lunar  Radio  Array  (LRA)  is  a  concept  for  a  radio  telescope 
operating  at  frequencies  around  100  MHz  and  sited  on  the  far 
side  of  the  Moon  with  a  prime  mission  of  making  precision  21- 
cm  cosmological  measurements.  With  reference  to  the  three  21- 
cm  cosmology  epochs  (3),  the  goal  for  the  LRA  is  to  probe  at 


3  http://www.skatelescope.org 
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least  the  epoch  of  first  star  formation,  and  as  deeply  into  the 
Dark  Ages  as  possible  (100  <  z  <  15,  10  MHz  <  v  <  100  MHz). 

Some  of  the  initial  designs  for  the  LRA  follow  the  config¬ 
uration  analysis  of  Lidz  et  al.  (2008),  who  advocate  a  highly 
compact  “super  core.”  For  such  a  configuration,  the  filling  fac¬ 
tor  is  large,  potentially  as  large  as  /  ~  0.9.  The  result  would  be 
to  reduce  the  required  integration  times  potentially  by  an  order 
of  magnitude.  However,  like  the  SKA-lo,  the  LRA  would  most 
likely  target  fields  well  away  from  the  Galactic  plane.  Similar 
comments  apply  to  the  LRA  as  to  the  SKA-lo  regarding  such 
observations. 

6.  Conclusions 

We  have  summarized  the  state  of  current  observations  and  the 
potential  for  future  observations  of  decame trie -wavelength  radio 
recombination  lines  (RRLs,  v  <  200  MHz).  These  lines  are  of 
interest  both  because  they  offer  a  sensitive  and  unique  probe  of 
cooler,  less  dense  ionized  regions  throughout  the  Galaxy  than  is 
probed  by  traditional  higher  frequency  RRL  studies,  as  well  as 
potentially  presenting  a  confusing  foreground  for  redshifted  21- 
cm  cosmology  and  astrophysical  measurements  from  the  Epoch 
of  Reionization  (EoR)  and  Dark  Ages. 

There  are  a  number  of  low  radio  frequency  arrays  ei¬ 
ther  under  construction  or  in  design  and  development,  includ¬ 
ing  the  Long  Wavelength  Array  (LWA),  the  LOw  Frequency 
ARray  (LOFAR),  the  low-frequency  component  of  the  Square 
Kilometre  Array  (SKA-lo),  and  the  Lunar  Radio  Array  (LRA). 
All  of  these  arrays  will  be  capable  of  observing  at  frequencies 
below  about  200  MHz,  and  they  could  be  used  to  study  RRLs, 
but  the  most  important  design  factor  for  a  telescope  is  the  areal 
filling  factor.  For  the  current  arrays,  LWA-1  has  a  filling  factor 
that  varies  with  frequency,  from  nearly  100%  at  its  lowest  fre¬ 
quencies  (»  20  MHz)  to  around  10%  at  the  highest  frequencies 
(~  80  MHz);  the  super-station  core  of  LOFAR  ranges  from  0.3 
to  0.02  over  a  similar  frequency  range  (its  LBA).  The  estimated 
integration  times  are  10  to  100  hr  to  detect  RRLs  at  optical 
depths  of  order  10  3;  in  the  case  of  the  LOFAR  super-station 
core,  these  relatively  modest  integration  times  are  achieved  in 
part  by  folding  potentially  hundreds  of  lines  over  large  band- 
widths  (up  to  48  MHz).  The  SKA-lo  and  the  LRA  are  still  in  the 
design  phase,  but  might  likely  have  higher  filling  factors,  poten¬ 
tially  larger  than  50%,  leading  to  the  detection  of  these  lines  in 
10  hr  or  less.  Further,  all  of  these  telescopes  offer  the  possibility 
of  improved  performance  as  a  function  of  frequency,  including 
both  improved  spectral  bandpasses  and  larger  frequency  tuning 
ranges. 

RRLs  have  been  studied  extensively  along  the  line  of  sight 
to  the  bright  radio  source  Cas  A.  However,  the  weakness  of 
the  lines  and  the  subsequent  need  for  long  integration  times 
and  careful  data  processing  to  remove  even  very  low  level  ra¬ 
dio  frequency  interference  (RFI)  have  hampered  efforts  for  the 
type  of  systematic,  multi-frequency  surveys  required  to  probe 
the  physics  of  the  cool,  low  density  ISM  throughout  the  Galaxy. 
Although  not  a  key  science  driver  for  any  of  these  telescopes, 
systematic  surveys  for  RRLs  could  be  important  secondary  sci¬ 
ence  to  emerge. 

From  the  standpoint  of  21 -cm  studies  of  the  Epoch  of 
Reionization  and  the  Dark  Ages,  deep  observations  (>  100  hr) 
will  likely  be  conducted  with  many  of  these  telescopes  to  search 
for  the  highly  redshifted  H  i  signatures,  implying  sensitivity  to 
much  lower  optical  depths  than  has  heretofore  been  possible, 
of  order  1CL4.  Simple  estimates  suggest  that,  if  there  is  RRL- 
hosting  gas  at  the  high  Galactic  latitudes  that  will  be  targeted  for 


the  observations,  the  RRLs  could  impose  spectral  fluctuations 
that  are  comparable  to  or  exceed  the  expected  H  i  signatures.  One 
significant  mitigating  factor,  however,  is  that  the  lines  are  likely 
to  be  quite  narrow,  relative  to  the  expected  Hi  features.  Hi  sig¬ 
natures  from  the  EoR  or  Dark  Ages  are  expected  to  have  widths 
of  order  0.2  MHz  or  larger;  by  constrast,  a  velocity  width  of  or¬ 
der  10  km  s  1  near  100  MHz  implies  a  line  width  of  only  a  few 
kHz.  Thus,  the  RRLs  might  likely  be  narrow,  with  a  low  spectral 
occupancy.  The  lines  might  easily  be  treated  as  non-terrestrial 
RFI. 

Two  final  comments  are  warranted  regarding  extremely  deep 
observations,  particularly  as  they  relate  to  efforts  to  detect  the  H  i 
signal  from  the  Epoch  of  Reionization  and  before.  It  is  not  yet 
known  what  is  the  most  relevant  frequency  range  for  21 -cm  cos¬ 
mology.  However,  it  is  possible  that  LOFAR-HBA  and  SKA-lo, 
and  possibly  the  future  LRA,  would  be  observing  in  the  regime 
for  which  RRLs  transition  from  emission  at  higher  frequencies 
to  absorption  at  lower  frequencies.  In  this  case,  the  lines  would 
not  be  detected.  Second,  most  of  our  comments  have  concerned 
^-transition  lines.  In  deep  integrations,  /?-,  y-,  and  potentially 
h-transition  lines  can  be  detected  (Stepkin  et  al.,  2007).  The  re¬ 
sulting  spectral  occupancy  would  be  much  higher  and  the  ability 
to  fold  or  remove  the  lines  could  be  reduced  substantially. 
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Table  1.  Known  Galactic  Carbon  a  RRL  Absorption  Lines" 


Sightline 

V 

(MHz) 

n 

T 

(xlO-3) 

^LSR 

(kms->) 

Av 

(km  s_1) 

Reference 

G0+0 

34.5 

575 

1.16 

5. 3+0. 8 

20.5+1.1 

5 

75 

443 

0.57+0.04 

-1.0+0. 5 

14.8+1.2 

APE88 

76 

441 

0.73 

-1 

24 

1 

GO-4* 

76 

441 

0.33 

-6 

13 

1 

GO-2 

76 

441 

0.57 

-10 

30 

1 

G0+0 

34.5 

575 

1.16 

5. 3+0. 8 

20.5+1.1 

5 

76 

441 

0.73 

-1 

24 

1 

GO+2 

76 

441 

0.70 

-2 

26 

1 

GO+4 

76 

441 

0.68 

-2 

31 

1 

G2-3.5* 

76 

441 

0.52 

7 

5 

1 

G2-2 

76 

441 

0.97 

5 

9 

1 

G2+0 

76 

441 

0.90 

2 

25 

1 

G2+2 

76 

441 

0.61 

6 

11 

1 

G3+0 

76 

441 

0.90 

-1 

14 

1 

G4+0 

76 

441 

0.54 

8 

17 

1 

G5+0 

34.5 

575 

0.74 

10.2+1.4 

21.2+2.0 

5 

G6+0 

76 

441 

0.73 

9 

25 

1 

G6.6-0.2 

76 

441 

0.87 

10 

28 

1 

G8+0 

76 

441 

1.09 

11 

22 

1 

G10+0 

34.5 

575 

0.81 

14.3+1.7 

37.5+2.5 

5 

0.93 

9.6+1.5 

24.9+2.3 

5 

0.47 

40.2+2.6 

18.4+3.7 

5 

76 

441 

0.72 

17 

26 

1 

G12+0.0 

76 

441 

0.91 

12 

20 

1 

G14-2.0 

76 

441 

0.88 

11 

32 

1 

G14+0 

34.5 

575 

0.66 

37.8+2.7 

54.0+3.8 

5 

76 

441 

0.85 

16 

25 

1 

G14+2.0 

76 

441 

0.51 

21 

28 

1 

G16+0.0 

76 

441 

0.76 

14 

47 

1 

M16 

68 

456 

1.8+0. 3 

18+2 

16+3 

6 

80 

435 

2. 0+0. 3 

20+1 

13+2 

6 

G16.5+0 

34.5 

575 

0.59 

26.4±2 

32.8+2.8 

5 

G06.9+0.8 

76 

441 

1.07 

20 

21 

1 

G18+0 

76 

441 

0.77 

15 

33+6 

1 

G20+0* 

76 

441 

0.59 

36 

18+3 

1 

G63+0 

34.5 

575 

0.42 

36.2+3.2 

46.0+4.5 

5 

G75+0 

25 

640 

1 

12 

15 

4 

G75+0 

34.5 

575 

0.40 

6.9+1. 8 

27.1+2.6 

5 

DR  21 

25 

640 

0.7+0. 3 

0+8 

42+12 

2 

34.5 

575 

0.74 

4.5+1.6 

18.8+2.2 

5 

S140 

25 

640 

1+0.3 

-6+17 

96+24 

2 

0.5+0. 3 

-36+17 

96+24 

2 

NGC  2024 

25 

640 

1 

10 

36 

4 

G287.4-0.6* 

76 

441 

0.48 

-23 

38 

1 

G312+0 

76 

441 

0.88 

-53 

16 

1 

0340+0* 

76 

441 

0.27 

-39 

30 

1 

G342+0* 

76 

441 

0.53 

-39 

17 

1 

G344+0 

76 

441 

0.65 

-27 

26 

1 

G346+0* 

76 

441 

0.20 

-20 

20 

1 

G348+0 

76 

441 

0.53 

-12 

14 

1 

G350+0 

76 

441 

0.81 

-11 

14 

1 

G352-2.0 

87 

441 

0.70 

-17 

23 

1 

G352+0 

34.5 

575 

0.65 

-1.2+1. 8 

34.8+2.6 

5 

76 

441 

1.25 

-10 

11 

1 
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Sightline 

y 

(MHz) 

n 

T 

(xlO'3) 

Llsr 
(km  s'1) 

Av 

(km  s'1) 

Reference 

G352+2.0 

76 

441 

0.83 

-8 

22 

1 

G354+0 

76 

441 

0.87 

-10 

26 

1 

G356+0 

76 

441 

0.76 

-4 

17 

1 

G358-2 

76 

441 

0.53 

-10 

28 

1 

G358+0 

76 

441 

0.90 

-3 

24 

1 

G358+2 

76 

441 

0.80 

3 

8 

1 

LI  407 

25 

640 

0.7 

-10 

17 

3 

(a>  This  listing  excludes  lines  for  the  Cas  A  line  of  sight. (hl  Tentative  detection 


References.  (1)  Erickson  et  al.  (1995);  (2)  Golynkin  &  Konovalenko  (1991a);  (3)  Golynkin  &  Konovalenko  (1991b);  (4)  Konovalenko  (1984); 
(5)  Kantharia  &  Anantharamaiah  (2001);  (6)  Anantharamaiah  et  al.  (1988). 
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Table  2.  Previous  and  Existing  Instrument  Capabilities 


Name 

Frequency 

(MHz) 

Resolution 

filling  factor 

* detect 

(hr) 

UTR-2 

26 

40' 

1 

i 

Gauribidanur 

34.5 

21'  x  25° 

0.59 

1.5 

DKR-1000 

40 

44'  x  1° 

0.2 

25 

Arecibo 

47 

80' 

0.25 

23 

Parke s 

76 

4° 

0.25 

2 

Notes.  rdeKct  is  the  time  to  detect  (5cr)  a  line  of  optical  depth  1CT3  with 
an  well-matched  instrumental  bandwidth.  The  instruments  listed  are 
ones  from  which  published  RRL  results  have  been  obtained.  In  many 
cases,  characteristics  have  been  taken  from  the  literature,  and  the  in¬ 
strument  itself  may  not  be  currently  operational. 


Table  3.  Future  Telescope  Capabilities 


Name 

Frequency 

(MHz) 

Resolution 

filling  factor 

T detect 

(hr) 

LWA-1 

20-80 

9°-2° 

0.9-0. 1 

-  10 

LOFAR-LBA 

30-90 

2°-0.6° 

~  0.1 

-  10 

LOFAR-HBA 

110-240 

30'— 15' 

~  0.1 

SKA-lo 

150 

7' 

~  0.1 

<  10 

LRA 

100 

10' 

~0.9 

<  10 

Notes.  Values  listed  are  characteristic  or  indicative,  but  none  of  these 
telescopes  have  entered  an  operational  state,  so  actual  performance  may 
vary.  For  LOFAR-LBA,  bandwidths  of  48  MHz  have  been  assumed. 
LOFAR-HBA  covers  a  poorly  explored  frequency  range  in  which  the 
lines  may  become  undetectable  as  they  transition  from  emission  at 
higher  frequencies  to  absorption  at  lower  frequencies;  the  same  may 
also  be  true  of  SKA-lo.  Particularly  for  the  SKA-lo  and  LRA,  final  de¬ 
sign  parameters  are  yet  to  be  determined,  and  these  values  are  indica¬ 
tive. 


Table  4.  Estimates  for  Detecting  Galactic  Plane  RRLs  with  LWA-1 


V 

(MHz) 

n 

/ 

^lines 

DF 

Ts  yS 

(xlO4  K) 

Afint 

(hr) 

20 

683 

0.92 

10 

0.26 

29.5 

50 

25 

640 

0.87 

8 

0.41 

18.3 

28 

30 

603 

0.76 

12 

0.59 

12.4 

6 

35 

572 

0.63 

10 

0.8 

8.9 

5.5 

40 

548 

0.5 

8 

1 

6.7 

7 

45 

526 

0.4 

6 

1 

5.2 

15 

50 

508 

0.33 

6 

1 

4.1 

23 

55 

492 

0.27 

10 

1 

3.3 

10 

60 

478 

0.23 

10 

1 

2.7 

14 

65 

466 

0.2 

8 

1 

2.3 

24 

70 

454 

0.17 

8 

1 

1.9 

32 

75 

444 

0.15 

6 

1 

1.7 

56 

80 

434 

0.13 

6 

1 

1.4 

72 

Notes.  Tsys  is  derived  from  the  GSM  sky  models  (de  Oliveira-Costa  et 
al.,  2008).  It  is  the  average  temperature  along  the  Galactic  Plane  towards 
Galactic  center  (-5°  <  b  <  5°,  340°  <  t  <  20°)  after  convolving  the 
models  to  LWA-1  beam  size  and  removing  the  pixels  with  the  highest 
and  lowest  2%  of  intensity.  The  dilution  factor,  DF,  assumes  4°  cloud 
sizes.  The  filling  factor  f  is  derived  using  the  station  dipole  layout  de¬ 
sign  (Kogan  &  Cohen,  2009)  and  accounts  for  overlaps  and  edge  effects. 
Mines  is  based  on  two  contiguous  bands  with  roughly  2  km  s-1  velocity 
resolution.  The  integration  time  to  detect  lines  with  optical  depths  of 
10~4  is  lOOx  the  values  listed  here. 
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